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ABSTRACT 

Context. HD 50138 is a southern star that presents the B[e] phenomenon, but its evolutionary stage is still not well known. This object 
presents spectral variability, which can be explained by outbursts or shell phases. Spectropolarimetric observations have shown the 
presence of a non-spherically symmetric circumstellar environment that is responsible for the B[e] phenomenon. However, up to now, 
the structure of the circumstellar medium of this object has not been studied deeply. 

Aims. Based on recent optical long baseline interferometric observations from the VLTI/MIDI and VLTI/ AMBER, and also from the 
Keck segment-tilting experiment, we study the structure of the circumstellar environment of HD 50138. 

Methods. The analysis of our data is based on geometrical analytical modeling, also using the recent LITpro software and considering 
a large space of parameters, which allows us to obtain a good estimate of the geometry of the circumstellar medium of HD 50138, 
responsible for the emission in the H, K, and N-bands. 

Results. We resolve the circumstellar environment of HD 50138 and describe its geometry for the first time in detail. Through anal- 
ysis of multiwavelength data, the presence of a dusty circumstellar disk with an orientation onto the sky-plane of 71 ±7°, which is 
perpendicular to the polarimetric measurements from the literature, was derived. We also derived that HD 50138 is seen under an 
intermediate angle related to the line of sight, 56 ± 4°. In addition, the structure of the disk and the flux contributions of the gas and 
dust components is discussed. 

Conclusions. Based on analysis of different sets of interferometric data, we describe the circumstellar disk whose geometric pa- 
rameters were determined, allowing us to understand the geometry of the circumstellar material of this bright star with the B[e] 
phenomenon. 

Key words. Stars: winds, outflows - Stars: individual: HD 50138 - Infrared: stars - circumstellar matter 



■ 1. Introduction 

The B[e] phenomenon in stars with quite different evolutionary 
stages (Lamers et al. '1998') is still a puzzle. Several questions, 
especially those related to the physical processes responsible for 
this phenomenon and associated to the role of binarity or even 
of stellar mergers, are still being debated. This happens because 
only a few objects have been studied deeply up to now and more 
than 50% of the B[e] stars have badly known or unknown evolu- 
tionary stages. 

In this paper, we continue the study of the southern galac- 
tic B[e] star HD 50138 (V743 Mon, MWC158, IRAS 06491- 



* Based on observations made with the ESO Very Large Telescope 
Interferometer at Paranal Observatory (Chile) under the Guarantied 
Time Observing, programs 080.C.0089, 080.D.0181 and 082.D.0047, 
and from the ESO archive. 

** This research made use of the Jean-Marie Mariotti Center LITpro 
service co-developped by CRAL, LAOG and FIZEAU. 



0654), started by Borges Fernandes et al. (120091 hereafter Paper 
I). This star has been studied since the beginning of the last 
century (Merrill et al. 1925 ), and several papers have described 
the presence of a strong spectral variability. This variability has 
been explained by the presence of outbursts and shell phases 
(Hutsemekers 119851 Andrillat & Houziaux I1991I I. Paper I has 
suggested that a new shell phase took place before 2007, based 
on the analysis of high-resolution spectroscopic data. The stel- 
lar parameters of this object were also derived, confirming that 
HD 50138 is a B6-7 III-V star 

Despite numerous papers on the topic, the evolutionary stage 
of HD 50138 remains highly uncertain. It has usually been clas- 
sified either as a pre main-sequence Herbig Ae/Be object or as 
a classical Be star (Paper I and references therein). This star has 
also been classified as "a transition object between a classical 
Be and a B[e] star" (Jaschek et al. ll993l l or simply as an unclas- 
sified B[e] star (Lamers et al. |1998l ). In addition, the possibility 
of binarity cannot be discarded (Cidale et al. 1200 II Baines et al. 
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I2OO6I 1. similar to what is found by Millour et al. (I2009I I for the 
unclassified B[e] star HD 87643. A deeper discussion about the 
nature of HD 50138 can be seen in Sect. 4 of Paper I. 

In Paper I, the color excess of HD 50138 due to purely in- 
terstellar extinction, i.e. E(B - V) = 0.08 mag, was also de- 
rived based on the existence of diffuse interstellar bands (DIBs). 
Moreover, from the analysis of the observed color indices, it 
was noted that the total extinction in the direction of this star 
was higher than the purely interstellar one. Therefore, this ex- 
tinction excess must come from the circumstellar environment 
of HD 50138. However, how gas and dust are distributed in the 
close vicinity of HD 50138 is still unknown. 

Measurements from polarimetry and spectropolarimetry 
have shown that the circumstellar material of HD 50138 is not 
spherically distributed. Based on UV and optical spectropo- 
larimetry and on the analysis of Ha intensity, Bjorkman et al. 
( 1199 8') show that electron scattering in a geometrically thin disk 
seen almost edge-on would be the main mechanism to produce 
the intrinsic polarization detected. On the other hand, Bjorkman 
et al. d 19981 1 also suggest that the dust responsible for the strong 
IR excess of HD 50138 (Allen [T973T l would be either distributed 
in an optically thin spherical envelope or formed by large grains 
(> lyum), not contributing to the intrinsic polarization measured. 
Paper I showed that an ionized circumstellar disk or a dusty op- 
tically thin envelope, composed mainly of small grains, could be 
responsible for the circumstellar extinction. However other sce- 
narios, like a dusty disk, where the dust is not located in the line 
of sight, could not be discarded. 

Bjorkman et al. (119981 1 find a polarization angle of 158°, 
which is quite similar to the mean value of 160+2.6° found by 
Yudin & Evans (T99E\ from UBV polarimetric measurements 
and of 157.5+5.0° by Oudmaijer & Drew (1999) from spec- 
tropolarimetric observations across the Ha region. In addition, 
Oudmaijer & Drew Cl999i) also suggest there are two different 
kinematic regions, where Ha would be formed: a rotating disk 
and also a more spatially extended component. Assuming a disk 
geometry, this polarization angle can be perpendicular or parallel 
to the major axis of the disk, depending on whether the circum- 
stellar material is optically thin or optically thick, respectively 
(AngelEMl Brown & McLean [19771 ). 

Recently, the advent of the optical long baseline interfer- 
ometry using large telescopes, especially the VLT 8-m and the 
Keck 10-m telescopes, has made it possible to resolve, and 
consequently describe, the gaseous and dusty circumstellar en- 
vironment of the brightest stars with the B[e] phenomenon 
(Domiciano de Souza et al. 2007; Monnier et al. 120091 Millour 
et al. 12009; Meilland et al. 2010). The high spatial resolution 
acquired, associated to the spectral dispersion provided by in- 
struments like VLTI/MIDI and VLTI/AMBER makes it possible 
to obtain sizes, shapes, and orientations of the circumstellar ma- 
terial, as a function of wavelength in the near and mid-IR ranges. 

In this paper, our study is therefore focused on the descrip- 
tion of the structure of the circumstellar material of HD 50138 by 
the analysis of recent interferometric data from VLTI/MIDI and 
VLTI/AMBER, and also from Keck segment-tilting experiment. 

The paper has the following structure. In Sect. 2 we present 
the observations and the data reduction processes. In Sect. 3, we 
present the reduced data and our analysis that considers geomet- 
rical models. In Sect. 4, we present the discussion of our results 
and in Sect. 5, a short conclusion. 



2. Observations and data reduction 

We obtained interferometric data on HD 50138 using the VLTI 
MID-infrared Interferometric instrument (VLTI/MIDI, Leinert 
et al. 2004). The observations were carried out between 
December 2007 and March 2009 (guarantied time observing) us- 
ing the 1.8-m auxiliary telescopes (ATs). The data were taken 
using the HIGH-SENS mode, where the photometry and the 
fringes are not recorded simultaneously. In total, 14 visibil- 
ity measurements were obtained considering different projected 
baselines and orientations onto the sky-plane, as described in 
the Table 1. Sirius and Aldebaran were used as calibrators in 
these observing runs and their angular diameters come from 
MIDI Calibrator Catalogued The (u,v) plane coverage of the 
VLTI/MIDI data of HD 50138 is seen in Fig. □ (left panel). 
The data reduction was performed using two different pack- 
ages: MIDI Interactive Analysis (MIA, Leinert et al. 120041 ) and 
Expert Work-Station (EWS, Jaffe l20D4b . version 1.5.1. Since 
the two packages provide similar results within the error bars, 
we decided to use the EWS reduced visibilities in this work. 
VLTI/MIDI allows us to obtain not only spectrally dispersed 
fringes, but also the low-resolution spectrum (R = 30) of the 
source in the N-band region (7.5-13 fim). A comparison be- 
tween VLTI/MIDI spectrum and IRAS low-resolution spectrum 
of HD 50138 is shown inFig.|2] 

In addition to our VLTI/MIDI data, we also used 92 N-band 
(at 10.7 jum) visibility measurements obtained from the Keck 
segment-tilting experiment It used 36 segments of 1.8-m each 
from the 10-m Keck telescope, forming four interferometric ar- 
rays, five subarrays, and baselines from 0.5-m to 9.0-m, by the 
tilting of the segments in the primary mirror (for more details, 
see Monnier et al. ,2009i) . The indication of the region covered 
by Keck data in the (u,v) plane can be seen in the Fig. [1] (left 
panel). 

We also analyzed sets of data at low spectral resolution 
(LR, i.e. R - 35) from five different observing nights in 
December 2007 and January 2008 from the Astronomical Multi 
BEam Recombiner, the near-infrared instrument of the VLTI 
(VLTI/AMBER, Petrov et al. l2003l) . We obtained data for two 
of these nights, and for the other three the data were taken from 
the ESO archive. A total of nine triplets, consequently measure- 
ments from 27 baselines, were available. However, based on the 
selection of the best quality data, we used only measurements 
from 17 baseline^ (Table 1). The (u,v) plane coverage of the 
VLTI/AMBER data of HD 50138 is presented in the Fig.[I](right 
panel). Measurements related to closure phases cannot be ana- 
lyzed due to the large error bars presented (about 50°). The cal- 
ibration stars were HD 52938, HD52436, HD96538, HD34137, 
HD44621, HD54810, and HD96068 (see Table 1), and their an- 
gular diameters come from Merand et al. (2005 ). The data re- 
duction was performed with the standard VLTI/AMBER data re- 
duction software (amdlib version 2. 1, Millour et al. 2004; Tatulli 
et al. i2007) and the frame selection was done as described by 
Cruzalebes et al. ( |2007] |. 

3. Geometrical models 

The analysis of our interferometric data is based on geometrical 
analytical modeling, which has allowed us to obtain a descrip- 

' http://www.eso.org/sci/facilities/paranal/instruments/midi/tools/ 
- Kindly provided by John Monnier 

^ The triplet taken on 2007-12-28, five baselines taken on 2008-01-06 
and two on 2008-01-09 were not used. Thus only VLTI/AMBER data 
from four different nights were actually used by us. 
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Table 1. VLTI/AMBER and VLTI/MIDI observing logs. 











projected baseline 






Number 


Date 


Obs.Time 


Stations 


Length 


PA 


Calibrators 


Ang. Diam. of the Calibrators 






(UT) 




[meters] 


[degrees] 


(used in the triplet) 


(mas) 


AMBER (ATs) 


Bi 


2007-12-29 


06:40 


Gl-DO 


62.3 


-36.4 HD 52436, HD 52938, HD 96538 


1.24, 0.90, 0.92 


B2 






DO-HO 


60.6 


74.2 






B3 






Gl-HO 


70.1 


17.7 






B4 


2008-01-02 


07:08 


KO-Gl 


90.3 


35.2 HD 34137, HD 44621, HD 52938 


0.83, 0.92, 0.90 


B5 






Gl-AO 


70.3 


-51.9 






Be 






KO-AO 


111.6 


74.2 






B7 


2008-01-06 


07:36 


Gl-AO 


62.3 


-44.5 


HD 52938 


0.90 


Bs 


2008-01-09 


05:26 


HO-GO 


31.4 


73.8 


HD 52938, HD 54810 


0.90, 1.49 


B, 






GO-EO 


15.7 


73.7 






Bio 


2008-01-09 


06:01 


HO-GO 


30.2 


74.2 


HD 52938, HD 54810 


0.90, 1.49 


Bii 






HO-EO 


45.3 


74.2 






B12 


2008-01-09 


06:43 


HO-GO 


27.9 


74.2 


HD 52938, HD 54810 


0.90, 1.49 


Bl3 






GO-EO 


14.0 


74.2 






Bl4 






HO-EO 


41.9 


74.2 






B,5 


2008-01-09 


03:08 


HO-GO 


29.7 


68.3 


HD 52938 


0.90 


B16 






GO-EO 


14.9 


68.3 






Bl7 






HO-EO 


44.6 


68.3 






MIDI (ATs) 


B, 


2007-12-09 


03:19 


Gl-HO 


68.9 


-12.3 


HD 48915 


5.85 


B2 


2007-12-09 


03:40 


Gl-HO 


68.6 


-10.2 


HD 48915 


5.85 


B3 


2007-12-12 


05:09 


Gl-DO 


71.6 


130.0 


HD 48915 


5.85 


B4 


2007-12-12 


06:00 


Gl-HO 


68.1 


6.8 


HD 48915 


5.85 


B5 


2007-12-13 


03:10 


Gl-DO 


66.6 


129.9 


HD 48915, HD 29139 


5.85, 19.90 


Be 


2007-12-26 


07:45 


Gl-DO 


58.1 


151.3 


HD48915 


5.85 


B7 


2008-11-09 


05:48 


EO-GO 


12.3 


60.7 


HD 48915 


5.85 


Bs 


2008-11-10 


05:47 


EO-HO 


37.1 


60.9 


HD 48915 


5.85 


B, 


2008-12-27 


06:35 


EO-GO 


15.5 


74.0 


HD 48915 


5.85 


Bio 


2008-12-28 


01:39 


EO-GO 


10.0 


50.8 


HD 48915 


5.85 


Bii 


2008-12-28 


02:36 


EO-HO 


36.9 


60.7 


HD 48915 


5.85 


B12 


2008-12-30 


01:35 


GO-HO 


20.3 


51.7 


HD 48915 


5.85 


B13 


2009-01-21 


06:08 


GO-HO 


27.1 


74.0 


HD 48915 


5.85 


Bl4 


2009-03-08 


01:22 


EO-HO 


47.7 


73.4 


HD48915 


5.85 



100 



50 - 



-50 - 



100 
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Fig. 1. The (u,v) plane coverage of HD 50138 obtained with VLTI/MIDI (left panel) and VLTI/AMBER (right panel). The region 
covered by Keck data is also indicated as inside of the stripped circle in the left panel. The dotted line in both panels represents the 
polarization angle measured by Yudin & Evans (I1998I I. 



tion of the spatial distribution of the flux from the source in dif- tained from Keck and VLTI/MIDI, which are mainly caused by 
ferent bands (H, K, and N-bands), considering a large space of circumstellar dust emission, 
parameters. The first step was to analyse our N-band data ob- 
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Table 2. VLTI/MIDI visibilities and equivalent Gaussian FWHM. 







Visibility 




Gaussian equivalent FWHM (mas) 


n° 


8yum 


10.7yum 


13yum 


8yum 


10.7//m 


13yum 


Baselines roughly parallel to the Yudin & Evans (2001) polarization measurement 


Bi 


0.10+0.06 


0.16+0.07 


0.30+0.12 


19.1+ 2.7 


22.4+ 2.8 


22.6+ 4.1 


B, 


0.11+0.06 


0.17+0.07 


0.35+0.09 


19.1+2.5 


22.4+ 2.6 


21.1+ 2.7 


B3 


0.07±0.05 


0.09+0.05 


0.22+0.06 


19.7+ 3.4 


25.1+ 3.4 


24.2+ 2.4 


B4 


0.04+0.04 


0.09+0.04 


0.42+0.10 


22.6+ 4.7 


25.9+ 2.4 


19.5+ 2.7 


B5 


0.09+0.06 


0.10+0.06 


0.21+0.09 


20.4+ 3.6 


26.0+ 3.3 


26.6+ 3.8 


Bs 


0.16+0.03 


0.30+0.04 


0.37+0.12 


20.2+ 1.0 


21.7+ 1.3 


24.4+ 3.9 



Baselines roughly perpendicular to the Yudin & Evans (2001) polarization measurement 



B7 


0.61+0.05 


0.62+0.05 


0.78+0.16 


50.3+4.3 


64.8+ 


5.7 


57.4+ 26.3 


Bs 


0.11+0.05 


0.14+0.06 


0.35+0.12 


34.7+ 3.7 


43.1+ 


4.5 


39.3+ 6.6 


B, 


0.54+0.08 


0.52+0.08 


0.88+0.13 


44.2+ 5.4 


60.1+ 


6.9 


33.5+ 24.9 


Bio 


0.51+0.08 


0.61+0.07 


0.89+0.23 


71.8+ 8.0 


80.3+ 


9.7 


49.4+ 46.4 


Bn 


0.12+0.07 


0.16+0.07 


0.44+0.13 


34.2+5.1 


42.1+ 


5.6 


34.8+ 6.4 


B [0 


0.53+0.10 


0.56+0.10 


0.75+0.16 


34.3+ 5.2 


43.2+ 6.5 


37.7+ 14.5 


Bl3 


0.30+0.04 


0.37+0.06 


0.55+0.08 


35.4+ 2.2 


41.9+ 


3.2 


40.2+ 5.1 


Bl4 


0.05+0.05 


0.07+0.05 


0.21+0.07 


32.1+ >11.8 


38.8+ 


6.1 


37.4+ 4.2 



80 F 




A (/.im) 

Fig. 2. VLTI/MIDI and IRAS LRS spectra of HD 50138. The 
wavelengths of crystalline silicate grains in the VLTI/MIDI 
spectrum range, e.g. ortho-enstatite, are overplotted as vertical 
solid Unes. 

3.1. N-banddata (Keck+VLTI/MIDI) 

3.1 .1 . The circumstellar environment at 1 0.7jum 

Monnier et al. (I2009I I obtained data at 10.7 yum with the Keck 
segment-tilting experiment by considering different orientations 
and baseline lengths from 0.5 to 9.0 m, making it possible to 
obtain information about the large-scale circumstellar environ- 
ment. They could model these data by assuming both (i) a 1-D 
Gaussian disk (actually a round 2-D geometry) with a full width 
half maximum (FWHM) of 58+6 mas, associated with an ex- 
tended circumstellar envelope (called "halo" by them) on scales 
larger than 0.5", which would be responsible for only 1% of 
the total flux in the N-band; and (ii) a 2-D Gaussian disk with 
FWHM for the major and minor axes of 66±4 and 46±9 mas, re- 
spectively, also in association with a "halo" responsible for just 
1% of the total flux in the N-band. The orientation of the ma- 
jor axis of the 2-D Gaussian disk was derived as 63+6°, which is 
roughly perpendicular to measurements of the polarization angle 
cited in the Sect. 1 . 



Thus, we have analyzed VLTI/MIDI data that have a broader 
spread of baselines in terms of lengths and orientations (Table 
1), allowing us to constrain both the radial profile of the object 
and its flattening. These data can be divided into two groups in 
terms of baseline orientations (Table 2): from Bi to Be, which 
are roughly parallel to the polarization angle, and from B7 to 
Bi4, which are roughly perpendicular to it. The parallel base- 
lines show an environment that is less resolved than seen by the 
perpendicular ones. This is shown in Table 2 from the estimate of 
the object extension (FWHM) assuming a simple 1-D Gaussian 
disk model. This result agrees with the orientation of the major 
axis of the disk found by Monnier et al. ( I2009I I. 

However, from our two groups of baselines, we derived a 
mean value for the FWHM of the major and minor axes at 10.7 
fim of 52+15 and 24±2mas, respectively. The large difference 
compared to the results of Monnier et al. (2009), i.e., a factor of 
1 .3 for the major and 1 .9 for the minor axis, shows the limitation 
of their calculations owing to the absence of longer baselines. 
In Fig.[3l we have plotted the VLTI/MIDI and Keck visibilities 
as a function of the spatial frequency, and we note that different 
single Gaussian disks with FWHM from 20 to 100 mas cannot 
model them. Since the trend in the visibilities is not compatible 
with other simple models, such as uniform disk or ring, we need 
to find a geometric analytical model that simultaneously fits the 
low spatial frequencies probed by the Keck and the higher ones 
measured by the VLTI/MIDI. We test two different models that 
are compatible with this quahtative analysis: 

1. A model consisting of two elliptical Gaussian distributions 
(representing a compact and a larger component) with the 
same flattening and orientation of the major axis (6). The 
other parameters are the FWHM of the major axis of each 
Gaussian and the relative flux of the compact component Fi 
(i.e., FWHM of the second Gaussian, F2 = 1 - Fi). Thus, 
this model has five free parameters. 

2. A flattened exponential distribution, i.e. I(r) cc s-\''/''o\. The 
model parameters are the flattening, the orientation of the 
major axis (6), and ro. To enable direct comparison with the 
Gaussian distribution, we use the FWHM of the exponen- 
tial distribution instead of iq. A simple calculation using the 
formula I(FWHM/2) = 0.5, direcdy gives FWHM = 1.38 ro. 
This model only has three free parameters. 



Borges Fernandes et al.: The galactic unclassified B[e] star HD 50138 



5 



Finally, we also use a single elliptical Gaussian distribution 
as a reference model. To not be biased by a putative change in the 
envelope geometry between 8yum and 13/im, we decided to first 
focus on the visibilities at 10.7//m. Using the VLTI/MIDI and 
Keck data, we have a set of 106 measurements. The model fitting 
is done by minimizing the reduced chi-square (x^.) on the visi- 
bilities, using a Levenberg-Marquardt (LM) two-dimensional al- 
gorithm. As this technique allows converging to local minimum, 
we use lO'* starting random positions in the parameter space to 
be sure that we reach the global minimum. The result of this 
model fitting is presented in Table |3] 

As already noticed in Fig. [3] a single elliptical Gaussian dis- 
tribution is not able to fit our dataset coiTectly, since its^^ is sig- 
nificantly higher than for the other two models. The 2-Gaussian 
model slightly shows the best x^., but the characteristic sizes, 
based on the FWHM of the major axis of the exponential distri- 
bution and of the compact component of the 2-Gaussian distri- 
butions, are compatible within 3cr, i.e., 30.7+1.4 and 35.2+1.5 
mas, respectively. Moreover, the flattening and orientation of 
the major-axes are consistent within Icr. In addition, we could 
also estimate the inclination related to the line of sight of the 
circumstellar structure, based on the flattening and assuming 
the 2-Gaussian and exponential distributions, i.e., 56.7+0.4° and 
59.0+1.6° respectively, which are also compatible within 2cr. 



Table 3. Parameters of the best-fit models from our 10.7yum anal- 
ysis. 



Gaussian 2-Gaussians Exponential 











0.68 ± 


0.04 






FWHM,„„,„, 


(mas) 


64.7 : 


to.e 


35.2 4 


: 1.5 


30.7 ± 


1.4 


FWHM,„„,„2 


(mas) 






131.4 4 


: 11.2 






Flattening 


2.94 + 


0.10 


1.82 ± 


0.02 


1.94 ± 


0.09 


en 




59.1 : 


t 1.7 


65.9 ± 2.0 


64.1 ± 


2.3 


in 




70.1 : 


t 0.7 


56.7 ± 0.4 


59.0 ± 


1.6 


xi 




5.1 


1.9 


24 





3.1.2. Wavelength dependence through the N-Band 

The VLTI/MIDI data show an object that is partly resolved for 
all baselines with a slight increase in visibilities with the wave- 
length. This mainly stems from the wavelength dependence of 
the interferometer spatial resolution, i.e. (/>,„,„ oc AfB, which is 
a factor 1.6 higher at S/^m than at 13//m, explaining the ob- 
served increase in visibility with the wavelength. No indication 
of a more extended envelope in the silicate band region is seen, 
unlike what is expected, since the opacity is supposed to be sig- 
nificantly higher in this band than in the continuum. However, 
as shown in Fig.|2] the silicate band is strongly in emission and 
is very broad, covering almost all the 8- 13/im VLTI/MIDI spec- 
tral band. This could explain why no clear drop is seen in the 
visibilities at the center of the silicate band. 

To quantitatively constrain the geometry of the circumstel- 
lar environment through the N-band, we used the same fitting 
method described in the previous section, considering only the 
2-Gaussian model, which provides the best results. For each 
VLTI/MIDI spectral channel, we fit the visibility of the 14 base- 
lines. Thus, the previously defined model parameters are now 
functions of the wavelength. In Fig. 31 we plotted the value of 
the best-fitted parameters as a function of the wavelength for the 
2-Gaussian model. 



We calculated an average over all spectral channels, i.e., 
< and found a value of 1.5. Moreover, as expected 

from the global trend of the spectrally resolved visibilities, the 
geometry does not strongly depend on the wavelength. The 
FWHM of major axis of the compact component varies from 
30+2 mas at Efj.m to 37±3 mas around lOjum and return to 
32+3 mas around 12 jum. It remains unclear whether these small 
variations, on the order of 2cr, are significant or not. This is also 
the case for the flattening of the compact component with vari- 
ations between 1.7 and 2.0, and for the orientation of the major 
axis with variations between 60° and 80°, which are both on the 
order of 2cr. The relative flux of the compact component is stable 
through the whole N-band; i.e., Fi = 0.7+0.03. Finally, the size 
of the extended component is larger than 100 mas. This means 
that it is fully resolved even by the shorter baselines. 

By analyzing the results in the frame of a circumstellar dusty 
disk, we can constrain a few geometrical parameters of the disk. 
First, we can deduce that the inner radius of the dusty disk lies 
within the compact component, so that its extension is smaller 
or equal to 15 mas. Considering a distance of 500 pc, this corre- 
sponds to an inner radius of the dusty disk smaller than 7.5 AU. 
Second, we can directly derive the inclination angle of the disk 
in the case of the geometrically thin disk; i.e., i = arccos (bi/ai) = 
57±3°. However, the opening angle of the disk may not be neg- 
ligible, so that we can, as shown by Meilland et al. 62007 ). only 
put a lower limit on the inclination of the object, i > 54°. Finally, 
we can also confirm that the major axis of the disk is roughly 
perpendicular to the polarization measurements, 9 = 70±12°. 

3.2. H and K-band data (VLTI/AMBER) 

The analysis of data from this spectral region is not trivial, since 
the near-IR emission usually comes from different sources: the 
stellar radiation, the hot gas located closer to the star than the 
dust sublimation radius (in a dust-free region), and also the hot 
dust, mainly located at the inner radius of the dusty disk. 

We started our analysis of the near-infrared VLTI/AMBER 
measurements by modeling the K-band data. The visibilities are 
plotted as a function of the spatial frequency in Fig |5] (right 
panel), where we clearly see that the object is fully resolved by 
the longest baselines. However, a non-negligible part of the flux, 
in the range of 10-20%, is still unresolved, which probably rep- 
resents the stellar flux. 

Another interesting point is the absence of a sharp edge, such 
as an inner rim of the dusty disk, which would produce oscil- 
lations in the visibilities. This absence can be for three differ- 
ent causes: (a) the inner radius is very extended and thus over- 
resolved even by the smallest baselines, (b) its contribution to 
the total flux is small enough to reduce the oscillation amplitude 
to an undetectable level, and (c) there is no gap in emission be- 
tween the inner gaseous envelope/disk and the inner rim of the 
dusty disk. We note that a combination of these causes can also 
be possible. Whatever causes the absence of oscillations, it will 
strongly prevent us from disentangling the gas and dust emission 
in the K-band. 

Thus, to perform a model fitting of these data, we decided 
to test LITprcQ (Lyons Interferometric Tool prototype), a new 
software developed by the Jean-Marie Mariotti Center (JMMC). 
It uses a set of elementary geometrical functions and also the 
LM algorithm to minimize thex^-- It is important to cite that we 
did not use LITpro for the N-band analysis, because it is not yet 
able to perform an achromatic model fitting. 

* LITpro software available at |http : //www.jmmc . f r/litprol 
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B/X (cycles/arcsec) 

Fig. 3. VLTI/MIDI N-band spectrally resolved visibilities plotted as a function of the spatial frequency (B/A) for the 14 baselines. 
Colors indicate the baselines position angle with respect to the polarization measurement by Yudin & Evans (I1998I I. The color 
scale (i.e., blue, cyan, green, yellow, orange, and red) goes from blue for baselines parallel to the polarization angle to red for the 
perpendicular ones. The dots with error bars indicate the 10.7/im Keck measurements. The thin solid and dashed lines represent the 
visibilities from Gaussian disks with different FWHM. The thick red solid line represents the model fitting aligned to the PA of the 
major axes of the 2-Gaussian distribution and the blue one to the PA of the minor axes. 



We decided to start our model fitting by using the simplest 
model that could fit the data, i.e. an unresolved star plus an el- 
liptical Gaussian model. Such a model can help us to derive the 
characteristic extension of the circumstellar environment. This 
simple model (hereafter Model 1) has the following parameters: 
the stellar relative flux (F*), and the Gaussian relative flux (Fi), 
its FWHM of the minor axis (FWHM„„„o,.i ), the orientation of 
the major axis on the sky-plane (0), and the flattening. Since the 
flux is normalized, we have Fi = 1 - F*, and consequently. Model 
1 has four free parameters. 

The result of the model fitting is presented in Table |4] 
and with = 40.8, we see that Model 1 failed to reproduce 
the spatial-frequency visibility variations (Fig. |5] right panel). 
However, we note that the flattening and orientation of the model 
1, i.e. 1.7+0.3 and 66+9°, are compatible within Icr with those 
ones derived from the lO.Vyum modeling. In addition, we could 
also, based on the flattening, estimate the inclination related to 
the line of sight of the circumstellar structure, i - 54+8°, which 
also agrees with our N-band results. 

Since the gas and dust emission regions may have different 
characteristic sizes, we also tried to model the environment with 
a more complex geometry, consisting of an unresolved star plus 
two Gaussian distributions, with the same flattening and orien- 
tation onto the sky plane. This model (hereafter Model 2) has 
the following parameters: the stellar flux (F*), the first Gaussian 
flux (Fi) and its FWHM of the minor axis (FWHM„„„o,.i), the 
second Gaussian relative flux (F2), its FWHM of the minor axis 
(FWHM„,„,or2), and the Gaussians flattening and orientation (0). 
As described for Model 1, owing to the normalized flux. Model 



Table 4. Results of the K-band model-fitting using LITpro. 



K-band 



Parameter 


model 1 


model 2 




0.19 ±0.03 


0.12 ± 0.01 


Fi 


0.81 ±0.03 


0.61 ± 0.06 


F2 




0.27 ± 0.07 


flattening 


1.7 ± 0.3 


1.7 ±0.3 


en 


66 ±9 


77 ±2 


FWHM„,„oH (mas) 


4.4 ± 0.5 


3.0 ± 0.4 


FWHM„„„,r2 (mas) 




> 14.0 


HI 


54 ±8 


54 ±8 


xi 


40.8 


13.3 



2 has six free-parameters. We obtained - 13.3, which is sig- 
nificantly better than the one for Model 1, but not yet fully satis- 
fying. 

We also tried to use other intensity distributions, such as uni- 
form disks, rings, and exponentials, but we did not manage to 
obtain a better fit. Thus, we decided to keep Model 2 as the ref- 
erence for the K-band (see Tableland Fig. |5] right panel). 

The orientation of the envelope, its flattening, and its inclina- 
tion related to the line of sight, determined by Model 2 (1.7+0.3, 
77±2°, and 54±8°, respectively) are also compatible with our 
N-band analysis. The most compact Gaussian major axis has an 
extension of 5.1 + 1.1 mas, corresponding to ~ 2.6 AU, whereas 
the extended one cannot be fully constraint, but it is larger than 
24 mas, i.e. 12 AU. The relative flux of the compact component 
(61+6%) dominates the K-band emission. 
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Fig. 4. Best-fit model parameters and for the 2-Gaussian model considering only the VLTI/MIDI data. All the parameters are 
plotted as a function of the wavelength: flux ratio of the compact component (top left); FWHM of the major and minor axes of the 
compact component (top center); FWHM of the major-axis of the large component (top right); flattening (bottom left); orientation 
of the major-axis (bottom center); and^^ (bottom right). 




Fig. 5. Model fitting of the VLTI/ AMBER data, using the parameters obtained by LITpro: right panel for the K-band assuming both 
scenarios as described in the text for Model 1 (dotted lines) and 2 (solid lines), and left panel for the H-band assuming a scenario 
like Model 2. The different colors follow the same description as in Fig. [3] representing the orientation of the baselines. 



As stated before, it is really hard to conclude about the sizes 
of the gas and dust emission regions from this simple model- 
ing. We know from the VLTI/MIDI modeling that the inner rim 
of the dusty disk is smaller than 7.5 AU, but both the compact 
(HWHM ~ 1 .3 AU) and extended components (HWHM ~ 6 AU) 
are compatible with this. Thus, we have tried to substitute in our 
models one of the components by a ring. Unfortunately, we did 
not manage to obtain any better agreement with such a model 
than for Model 2. Nevertheless, the ring radius cannot be smaller 
than 3 AU, otherwise oscillations in the visibilities would be 
seen even for the longest baselines. Thus, it finally seems that 
the large component probably represents the inner radius of the 
dusty disk. 



In addition, we have also tried to model the H-band 
VLTI/ AMBER data, where the quality is clearly lower compared 
to the K-band data (Fig. |5] left panel). Using the same scenario 
as described by Model 2, we did not manage to constrain prop- 
erly the model parameters, especially the extensions, flattening, 
and orientation of the circumstellar environment. However, we 
can consider that only the relative flux of the components would 
change in the H-band compared to the K-band. Thus, we fixed 
the geometrical parameters with the same values as used for the 
K-band, and we derived the fluxes for the three components: F* 
= 0.22 + 0.04, Fi = 0.65 + 0.10, and F2 = 0.13 + 0.07. The 
corresponding visibilities are overplotted in Fig. |5] (left panel). 
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4. Discussion 

A good coherence was obtained considering the model-fitting 
of VLTI/ AMBER, VLTI/MIDI and Keck data, allowing us to re- 
solve and describe in detail for the first time the circumstellar en- 
vironment of HD 50138. In addition, based on our data and their 
different project baselines and position angles, we cannot ver- 
ify that there is variability in the interferometric measurements 
associated to the spectral variations of this object that were de- 
scribed in Paper I. Thus, we can now draw a general picture for 
this object. 

4.1. General structure of the dusty disk 

Based on our model fitting for different data, we saw that the 
star is seen under an intermediate inclination angle of 56+4° and 
is surrounded by a large dusty disk, with a major-axis roughly 
perpendicular to the polarization angle, 6 = 71+7°). This disk 
probably has a small opening angle, which can explain the pho- 
tospheric lines seen in the high-resolution spectra of HD 50138 
(Paper I). Moreover in this case, since the dust would be out of 
our line of sight, it does not significantly contribute to the visi- 
ble polarization, in agreement with the results of Bjorkman et al. 

From our models in the mid-IR, we see that the dusty disk 
emits up to a large distance from the central star, i.e. more than 
30 AU. Moreover, the flux measured by the VLTI/MIDI spec- 
trum is about 13% lower than the value provided by IRAS LRS 
(see Fig |2]l, both with similar spectral resolutions. This differ- 
ence can be considered quite small, since the uncertainty of the 
VLTI/MIDI flux is around 20%. In addition, the VLTI/MIDI 
field of view (FOV) is around 1" at 8/im and of 1.5" at 13/im, 
while the IRAS LRS one covered a region of 7.5'. Based on this, 
we can consider that VLTI/MIDI spectrum represents the whole 
N-band emission of HD 50138, thereby avoiding any external 
contamination. 

4.2. Mineralogy of HD 50138 dust 

The VLTI/MIDI spectrum also allows us to qualitatively analyze 
the mineralogy of the circumstellar medium of HD 50138. We 
note the presence of a lO/im emission bump that corresponds 
to the silicate band, which does not seem to affect the visibil- 
ity measurements. However, as already mentioned in Sect. 3.1.2, 
this is probably because the silicate band is as broad as the 
VLTI/MIDI spectral bandwidth. 

In addition, the VLTI/MIDI spectra, even from shorter base- 
lines, show features from cristalline silicate grains, such as 
ortho-enstatite (Fig. |2]i. This does not agree with the expecta- 
tions for young stellar objects (van Boeckel|2008]and di Folco et 
al. |2009l l. However, a more detailed analysis of the dust compo- 
sition, associated to a larger sample of data, becomes necessary 
to provide a complete study of the dust content, as a function of 
the distance from the star. 

4.3. Separating the stellar, gas, and dust emissions 

In our K-band modeling, we found that HD 50138 circumstel- 
lar environment can be divided into two components. The first 
one, fully resolved even by the smallest baselines, contributes to 
27+7% of the flux. A lower limit on its extension around 24 mas 
(12 AU) can be set. Comparing this component with the N-band 
modeling, we conclude that it probably represents the K-band 
contribution of dusty disk inner rim. However, only a radiative 



transfer code, which couples gas and dust, will allow us to derive 
the radius of the inner dusty disk properly. It is also clear that, 
even in this case, we will not find an specific inner radius, but a 
range of radii, depending of the grain size. 

The second component is much more compact with FWHM 
- 5.1 + 1.1 mas (i.e. 2.6 AU), and it represents 61+6% of the K 
band, which can be associated to hot gas emission. Concerning 
the mid-IR emission, since the N-band stellar flux itself is ob- 
viously negligible, it can only originate in the dust surrounding 
HD 50138. 

Through the stellar, dust, and circumstellar gas relative 
fluxes, and photometric and spectro-photometric measurements 
in H, K, and N bands, it is possible to derive the flux contribution 
of each component. We used 2MASS measurements to estimate 
the total flux for the H and K-bands, and the VLTI/MIDI spec- 
trum (as cited in the Sect. 4.1) for the total flux of the N-band. 
The resulting fluxes are shown in Table |5] 

4.4. The gaseous inner disk 

Thanks to the previous section, we now have suggested that the 
compact component of our K-band model represents the emis- 
sion of the hot gas surrounding HD 50138. Moreover, the flatten- 
ing of this component, i.e. 1.7+0.3, is compatible with the one 
determined in N-band for the dusty disk, i.e. 1.82+0.02, consid- 
ering the geometrical model at 10.7/vm. The major axis orien- 
tation of these two components is also compatible. Therefore, 
most of the hot gas must be located in a disk. This disk may be 
a Keplerian accretion disk, but just in case of a young nature for 
HD 50138, or even an outflowing disk wind for a more evolved 
scenario. 

To distinguish between these different scenarios, it is neces- 
sary to constrain not only the circumstellar environment geom- 
etry but also its kinematics. This can be done by coupling the 
high-angular resolution of modern interferometers with a suffi- 
ciently high spectral resolution. Based on this, observing pro- 
grams for VLTI/AMBER are currently being prepared, using its 
high spectral resolution mode (R=12000). 

Another interesting point is to compare the assumption of a 
geometrically thin disk seen under an intermediate inclination 
angle with the Balmer lines profiles described in Paper I, i.e. 
Ha, H/3, Hy, and H6. In spite of their variations, all these emis- 
sion lines clearly show a strong and narrow absorption compo- 
nent that is probably formed in the circumstellar environment. 
Thus, a non-negligible part of the emitting gas must be located 
in the line of sight, i.e 56°, hence far above the equatorial disk. 
Consequently, this gaseous emission could come from the stellar 
wind, as described for Herbig stars by Kraus et al. (I2008I I. 



Table 5. Fluxes of the components of HD 50138 for the different 
bands. 



Band 


A 




Fluxes in W m um ' 






(yum) 


Star 


Gas 


Dust 


Total 


H 


1.65 


2.4 10"'^ 


7.0 10^'^ 


1.4 lO^" 


1.1 10"" 


K 


2.20 


1.1 10^'^ 


5.6 10-'^ 


2.5 10^'^ 


9.2 10^'^ 


N 


8 






2.3 10-'^ 


2.3 10-'^ 


N 


10.7 






1.6 10-'2 


1.6 10-'2 


N 


13 






6.4 lO-'-'' 


6.4 lO-'-'' 
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5. Conclusions 

HD 50138 is a curious star whose evolutionary stage is un- 
known. In this paper, we could describe the geometry of both 
the gaseous and dusty circumstellar environments, based on the 
analysis of recent interferometric data. Through geometrical an- 
alytical modeling, it was possible to resolve a circumstellar disk 
around HD 50138 and to describe its parameters for the first time 
in detail. 

However, more observations, especially from the VLTI and 
with high S/N, are planned. The analysis of these data, in as- 
sociation with 3-D codes, will allow us to derive not only the 
physical parameters of the dusty environment precisely, but also 
those from the gas emitting regions. In addition, new data wiU 
make the image reconstruction of HD 50138 possible. This in- 
formation will be extremely valuable for better understanding 
the physical mechanisms responsible for the formation of the 
B[e] phenomenon. 
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